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Abstract 

Using a large sample of D + — ► K~ir + ii + v decays collected by the FOCUS photo- 
production experiment at Fermilab, we present the first measurements of the he- 
licity basis form factors free from the assumption of spectroscopic pole dominance. 
We also present the first information on the form factor that controls the s-wave 
interference discussed in a previous paper by the FOCUS collaboration. We find 
reasonable agreement with the usual assumption of spectroscopic pole dominance 
and measured form factor ratios. 



1 Introduction 

The D + — > K~it + n + v decay is described in terms of helicity basis form factors 
that give the q 2 dependent amplitudes for the K to be in any of its possible 
angular momentum states [1]. Traditionally [2,3], these helicity basis form 
factors are written as linear combinations of vector and axial form factors 
that, in turn, are assumed to have a q 2 dependence given by spectroscopic 
pole dominance. The pole masses are fixed to the known masses of the excited 
D+ states with vector and axial quantum numbers. This paper uses a new 
weighting technique to disentangle and directly measure the q 2 dependence of 
these helicity basis form factors free from the assumption of spectroscopic pole 
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dominance. We believe this paper represents the first non-parametric analysis 
of the D + — > K~TT + fi + u helicity basis form factors. 

Five kinematic variables that uniquely describe D + — > K~n + fi + u decay are 
illustrated in Figure 1. These are the K~ir + invariant mass (m^), the square 
of the [iv mass (g 2 ), and three decay angles: the angle between the tt and 
the D direction in the K~tt + rest frame (#v), the angle between the v and 
the D direction in the fiu rest frame (9e) , and the acoplanarity angle between 
the two decay planes (x)- The full intensity distribution for D + — > K~7i + fi + i / , 




Fig. 1. Definition of kinematic variables. 

differential in these five kinematic variables, is given in Reference [2]. Through- 
out this paper, we will use the simpler expression given by Eq. 1 that gives 
the form of this decay intensity after integration over the acoplanarity angle 
X- The x integration significantly simplifies the intensity by eliminating all 
interference terms between different helicity states of the virtual W + with rel- 
atively little loss in form factor information. 1 Following the notation of [1], 
the Eq. 1 intensity is written in terms of the four helicity basis form factors: 
H + (q 2 ), H (q 2 ), H_(q 2 ), H t (q 2 ) along with an additional form factor h(q 2 ) that 
describes the coupling of an additional, small s-wave amplitude contribution 
(with a constant amplitude Ae tS ) that is discussed in Reference [4]. Our ear- 
lier papers [4,2] assumed that the s-wave amplitude had the same zero helicity 
form factor as that for D + — > K*°fi + u - e.g. h (q 2 ) = H (q 2 ). 2 The acopla- 
narity averaged intensity in terms of the decay angles, helicity form factor 
products, and K Breit-Wigner amplitude (BW) is: 



((1 + cos0*) sin# v ) 2 |# + (g 2 )| 2 |BW| 2 
-((1 - cos0*) sin# v ) 2 |#-(g 2 )| 2 |BW| 



f t AM, q-mi +U1 -cos^jsmtfyj -"-W 
/ lA^dx = 1 { 

J 8 +(2sin^cosfl v ) 2 |#o(g 2 )| 2 |BW l 

^ +8 sin 2 d t cos6 v H (q 2 )h o (q 2 )Re{Ae- i5 BW} t 

1 The acoplanarity acceptance in our spectrometer is uniform to within 2%, hence 
no correction is required given the size of our statistical errors. 

2 Eq. 1 drops the term which is second order in the small amplitude A. 
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(1) 



IBW 



- ^-^{q 2 -n4)-^{ 



(sin^ sin# v ) 2 |#+| 2 + (sin0< sin 6> V ) 2 |^T_| 

+(2cos0< cosfl v ) 2 |# | 2 

+ (2 cos 0y) 2 1 H t | 2 + 8 cos fl £ cos 2 v #o#t 



where 



BW 



* \ 3 



m 2 Kn - m 2 , + im T (j^) 



(2) 



The first term gives the intensity for the /i + to be right-handed, while the 
(highly suppressed) second term gives the intensity for it to be left-handed. 3 

Before describing the non-parametric approach, we begin with a description of 
the traditional experimental approach. The D + — > K decay amplitude 

is typically analyzed [1] in terms of four form factors. This intensity expression 
is written in terms of four helicity basis form factors that are in turn written 
as linear combinations of vector and axial form factors as given in Eq. 3. 



H ± (q 2 ) = (M D + m K7r )A 1 (q 2 )T2 



M D K 



V(q 2 ) , 



Ho{q 2 ) 



HM 2 ) 



2m KnVq 2 



M D K 
rriKirVq 2 



{M 2 D - m 2 Kn - q 2 ){M D + rriK^A^q 2 ) 



-4- 



M 2 D K 2 
~M D + m K7T 



A 2 {q 2 



(3) 



(M D + m Kn )Mq 2 ) - mL + ?2 W ) 



M D + m K7T 



2q 2 



M D + m K7r 



A 3 (q 2 ) 



where K is the momentum of the K~n + system in the rest frame of the 
D + . The vector and axial form factors are generally parameterized by a pole 
dominance form: 



AM 2 ) = 



A(0) 



q 2 /M 



and V{q 2 



V(0) 



q 2 /M, 



v ' 



(4) 



3 We are using a p-wave Breit-Wigner form with a width proportional to the cube 
of the kaon momentum in the kaon-pion rest frame (P*) over the value of this 
momentum when the kaon-pion mass equals the resonant mass {Pq). The squared 
modulus of our Breit-Wigner form will have an effective P* 3 dependence in the 
numerator as well. Two powers P* come explicitly from the P* in the numerator of 
the amplitude and one power arises from the 4 body phase space. 
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where all previous experiments used spectroscopic pole masses of Ma = 2.5 GeV/c 2 
and My = 2.1 GeV/c 2 which in the case of D + — > K*°fi + u are tied to masses 
of vector and axial D* states. 

Using the spectroscopic pole dominance assumption, previous experiments [2,3] 
have fit the shape of the D + K intensity to at most 3 parameters 

which are ratios of form factors taken at q 2 = : r v = V(0)/Ai(0), r 2 = 
A 2 (0)/A 1 (0), r 3 = A 3 (0)/Ai(0) and in one case [2] the (constant) s-wave 
complex amplitude. 

As in our earlier paper [5] on the q 2 dependence of the D° —> K~fi + v form 
factor, we present the first non-parametric measurements of the helicity basis 
form factors that describe D + — > K^7i + fi + u decay. In particular, we will 
provide information on H±(q 2 ), H^q 2 ), and h (q 2 ) x H (q 2 ) in bins of q 2 by 
projecting out the associated angular factors given in Eq. 1. The cross term 
ho(q 2 ) x H (q 2 ) represents the interference between the s-wave and the K*°. 

Throughout this paper, unless explicitly stated otherwise, the charge conjugate 
is also implied when a decay mode of a specific charge is stated. 



2 Experimental and analysis details 

The data for this paper were collected in the Wideband photoproduction ex- 
periment FOCUS during the Fermilab 1996-1997 fixed-target run. In FOCUS, 
a forward multi-particle spectrometer is used to measure the interactions of 
high energy photons on a segmented BeO target. The FOCUS detector is a 
large aperture, fixed-target spectrometer with excellent vertexing and parti- 
cle identification. The FOCUS experiment and analysis techniques have been 
described previously [4,6,7]. 

To isolate the D + — » K~n + topology, we required that candidate muon, 
pion, and kaon tracks formed a secondary vertex with a confidence level ex- 
ceeding 25%. We required a primary vertex consisting of at least two charged 
tracks. The muon track, when extrapolated to the shielded muon arrays, was 
required to match muon hits with a confidence level exceeding 5%. The kaon 
was required to have a Cerenkov light pattern more consistent with that for a 
kaon than that for a pion by 1 unit of log likelihood [7]. No Cerenkov require- 
ment was made on the pion. 

To further reduce muon misidentification, a muon candidate was allowed to 
have at most one missing hit in the 6 planes comprising our inner muon sys- 
tem and a momentum exceeding 10 GeV/c. In order to suppress muons from 
pions and kaons decaying within our apparatus, we required that each muon 
candidate had a confidence level exceeding 2% to the hypothesis that it had 
a consistent trajectory through our two analysis magnets. 
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Non-charm and random combinatoric backgrounds were reduced by requiring 
both a detachment between the vertex containing the K^n + fi + and the pri- 
mary production vertex of at least 10 standard deviations and a minimum 
visible energy (E K + E w + E^) of 30 GeV. To suppress possible backgrounds 
from higher multiplicity charm decays, we isolate the Kit\i vertex from other 
tracks in the event (not including tracks in the primary vertex) by requiring 
that the maximum confidence level for another track to form a vertex with 
the candidate be less than 0.1%. 

In order to allow for the missing energy of the neutrino in this semileptonic 
D + decay we required the reconstructed Kir^x mass be less than the nominal 
D + mass. Background from D + — > K~tt + tt + , where a pion is misidentified 
as a muon, was reduced using a mass cut: we required that when the muon 
track is treated as a pion and the combination is reconstructed as a Kirn, 
the Kirir invariant mass was less than 1.820 GeV/c 2 . In order to suppress 
background from D* + — > D°ir + — > (K~ fj, + iy)n + , we required M(K~ /j, + i/ir + ) — 
M(K~fi + v) > 0.18 GeV/c 2 . The wrong-sign subtracted ttik-k distribution for 
these D + — > K~n + fi + u candidates is shown in Figure 2. Wrong-sign events 
have tracks identified as K~n + ji~ isolated in a detached vertex. 
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Fig. 2. Wrong-sign subtracted D + — > K~ir + ji + u signal. Over the full displayed 
mass range there is a right-sign excess of 14 798 events. For this analysis, we use 
a restricted mass range from 0.846-0.946 GeV/c 2 (shown by vertical lines). In this 
restricted region, there is a right-sign excess of 11 397 events. 



We will use a restricted mass range from 0.846-0.946 GeV/c 2 in this analysis 
in an effort to diminish the dependence of the helicity basis form factors on 
m Kn through Eq. 3. 

The technique used to reconstruct the neutrino momentum through the D + 
line-of-flight and tests of our ability to simulate the resolution on kinematic 
variables that rely on the neutrino momentum are described in Reference [4]. 
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3 Projection Weighting Technique 



In this section, we describe the weighting technique that we use to extract 
the helicity basis form factors that describe the D + — > K~ir + /j, + i/ angular 
distribution according to Eq. 1. For a given q 2 bin, a weight will be assigned 
to the event depending on its 8y and Op decay angle. We consider 25 joint 
A cos 8y x A cos 0£ angular bins: 5 evenly spaced bins in cos 6 V times 5 bins 
in cos6>£. Each event will acquire a weight designed to project out a given 
helicity form factor that depends on which of the 25 angular bins that it is 
reconstructed in. It is convenient to think of weighting as constructing a dot 
product of the form P a -D where D is a data vector that consists of the number 

of events reconstructed in each of the 25 angular vectors D = (m n 2 ... n 2 5 ) 

— * 

and P a is a projection vector for the a helicity form factor. The 25 components 
of each P a vector give the weights applied to each event reconstructed in one 
of the 25 angular bins. Eq. 5 says the product P + - D is equivalent to weighting 



the events in angular bin 1 by 



, weighting the events in angular bin 2 by 



etc. 



P+-D 



Ui + 



n 2 + 



25 



n 2 5 ■ 



(5) 



The P a weights are designed to project out the helicity form factors using 
Monte Carlo inputs. Here is how they are obtained. To simplify our discussion, 
consider the case of just three form factors H^_(q 2 ), H^q 2 ), H^q 2 ). For each 

q 2 bin, let {fh a } = ( fh + rh^ m ) where rh a is the number of events present in 

each of the 25 angular bins when a nominal H a (q 2 ) is turned on and all other 
nominal Ha^p^q 2 ) are turned off. As indicated in Eq. 6, for each q 2 bin the 

— * 

Di vector can be written as a linear combination of the three m vectors with 
coefficients f a {q^)- The f a {qf) functions are proportional to the true H 2 (q 2 ) 
along with pre-factors such as q 2 — m 2 and corrections such as acceptance and 
resolution. 



A = f+{q 2 i) rn + + f_{q 2 ) m_ + f (q 2 ) m . 



(6) 



We can convert Eq. 6 to the "component equation" shown in Eq. 7. 



m + - Di 

— * 

m_ • Di 
K rh - Di J 



^ m + ■ m + m + ■ m_ m + ■ m ^ 



(7) 
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The solution to Eq. 7 can be written as: 

U(q 2 ) = l P + ■ A , f-(qt) = l P~ ■ A , /ofe 2 ) = l Po ■ A , 



where our notation is that the left superscript gives the q 2 bin number and 
the subscript identifies the helicity of the projector. 



The l P a weights are given by Eq. 9. 



i p 



m + ■ m + m + ■ m_ m + ■ m 
m_ • m + m_ • m_ m_ • mo 



m + 1 



(9) 



One can correct for acceptance by using the proportionality relations such as 
Eq. 10, 



^ . ^ ^ (Hjq 2 )) 2 ■ - - 



(^ + (% 2 )) 2 



l H 4 



l P+ ■ Mi 



A = V+ • A , 



(10) 



where Mj are the bin populations from a Monte Carlo, generated assuming a 
trial form factor set H+(q 2 ), iJ^.(g 2 ), and H^q 2 ). As indicated in Eq. 10, the 

projection weights l P + and the projection-weighted Monte Carlo distributions 
( l P + ■ Mj) can then be used to construct an adjusted weight vector l p + . The 
(arbitrarily normalized) form factors H^(q 2 ), H 2 (q 2 ) and H^q 2 ) would then 
be obtained by making three weighted histograms using the t p + , % p_ and l p 
weights respectively. 

We next discuss some of the complications in applying the projective weighting 
scheme in our experiment due to smearing of kinematic variables because 
of the missing neutrino. In the absence of substantial q 2 smearing, a totally 
arbitrary set of trial form factors (such as H a (q 2 ) = 1) can be used to get 
unbiased estimates of the helicity basis form factors. This is no longer true 
when smearing in the kinematic variables q 2 , cos#v, cos^ is substantial since 
one must use reconstructed quantities in dealing with the data. The acceptance 
and resolution that affect the l p + , % p_ and *p weights depend on the true 
kinematic variables. Since the mapping between the true kinematic variables 
and reconstructed kinematic variables depends on the underlying form factors, 
one can bias the returned form factors to the extent that H a (q 2 ) ^ H a (q 2 ). 
We found through Monte Carlo simulation that it is only important to get a 
reasonably good first guess for H a (q 2 ) and multiple iteration is not required 
given the size of our statistical error bars. 
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Although the mass terms of Eq. 1 are suppressed by two powers of the muon 
mass, they are surprisingly important at low q 2 . Our analysis includes projec- 
tors for each of the six form factor products present in Eq. 1. Although we were 
unable to obtain useful information on H 2 {q 2 ) or the H (q 2 ) x H t (q 2 ) interfer- 
ence term, we needed to allow for them in the construction of l P a using Eq. 9 
to insure that the projectors used for H±(q 2 ), H^q 2 ), and h (q 2 ) x H (q 2 ) 
will be "orthogonal" to the angular terms associated with the H 2 (q 2 ) and 
Ho{q 2 ) x H t {q 2 ) contributions of Eq. 1. Without the incorporation of the mass 
term projectors, we see a dramatic mismatch between the input and output 
form factors for H±(q 2 ) in the first q 2 bin in our Monte Carlo studies. For 
example, if^(g 2 ) nearly drops to zero in the first q 2 bin if projectors for the 
mass terms are not included. 



Figure 3 summarizes a complete Monte Carlo simulation of the projective 
weighting technique. This Monte Carlo was run with 9 times our data sample 
but we have inflated the error bars by a factor of three to indicate the estimate 
of the errors expected in the data. The form factor measurements are plotted 
at the abscissa of the average generated q 2 for each of the 6 evenly spaced 
measured q 2 bins rather than at the measured bin center. The good agreement 
between the input and output form factors validates our assumption that the 
\A\ 2 term can be dropped when constructing the projective weights since these 
terms are included in the Monte Carlo simulation. 




).0 0.2 0.4 0.6 0.8 1 



0.2 0.4 0.6 0.8 1.0 



q 2 GeV/c 2 



Fig. 3. Monte Carlo study of the projective weighting technique. The form factors 
used in this Monte Carlo simulation are shown as the solid curves. The reconstructed 
form factors are the points with error bars. They are plotted in "arbitrary" units 
but the same unit is used for all four form factors in order to convey the relative 
size of each contribution. The plots are: (a) H^(q 2 ), (b) H^^q 2 ), (c) Hq^ 2 ), and 
(d) h (q 2 ) x H (q 2 ). 



The method does an excellent job at reproducing the input form factor prod- 
ucts (shown as a curve) both in terms of the shape and relative contribution 
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from each of the 4 form factor products. The input form factors used the 
measurements and s-wave model developed in Reference [2]. 

Figure 4 and Table 1 show the results obtained for data. The form factor mea- 
surements are plotted at the abscissa of the average generated q 2 for each of 
the 6 evenly spaced measured q 2 bins as determined from the Monte Carlo. 
The horizontal error bars given the r.m.s q 2 resolution for each bin. We have 
subtracted the weighted distribution obtained for background events simu- 
lated in our charm Monte Carlo that incorporates all known charm decays. 
Non-charm backgrounds are primarily eliminated through a wrong-sign sub- 
traction. Applying more stringent cuts in the analysis indicates that the only 
significant deviation from the model occurs in plot a) of Fig. 4 in the lowest 
bin. While we believe the low bin in H + (q 2 ) represents a very small fraction 
of the cross section for the decay and is, hence, particularly susceptible to 
unanticipated backgrounds, we cannot rule out that this effect originates in 
the physics of the real decay. 

We will refer to this representation of the data as the un- convolved analysis 
since the vertical error bars in Figure 4 do not represent the uncertainty in the 
form factors averaged within each bin boundary. This is because we have not 
corrected for the considerable smearing between the various q 2 bins by using 
the deconvolution technique [5] discussed in the next section. To the extent 
that the form factors vary smoothly, the un-convolved representation should 
still be faithful to the underlying form factors as was the case of the Monte 
Carlo study shown in Figure 3. We include the un-convolved results since they 
show more q 2 bins than possible in our deconvolution analysis discussed in the 
next section. 

0006 
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0.0000 u .™ 

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 
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0.0025 
0.0000 

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

q 2 GeV/c 2 

Fig. 4. Non-parametric form factors obtained for data with horizontal error bars 
given by our r.m.s. resolution. The plots are: (a) H^(q 2 ), (b) H 2 _{q 2 ), (c) H^q 2 ), 
and (d) /io(<? 2 ) x H^q 2 ). The r.m.s q 2 resolution varies (monotonically) from 0.125 
GeV 2 /c 4 , in the first bin, to 0.218 GeV 2 /c 4 in the last bin. 
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Table 1 

Summary of un-convolved results. This is a tabular summary of the data of Fig. 4 
multiplied by 1000. 



a 2 


a 2 


H 2 


H 2 


Hi 


ICQ A 11Q 


0.169 


r~i -I OCT 

0.125 


0.42 ± 0.09 


2.02 ± 0.87 


8.91 ± 


1.71 


4.3b ± 1.53 


0.272 


0.136 


0.31 ± 0.07 


1.75 ± 0.22 


4.09 ± 


0.54 


f 1 r\ to 

2.53 ± 0.78 


0.400 


0.154 


0.32 ± 0.05 


1.62 ± 0.14 


2.44 ± 


0.32 


1.03 ± 0.57 


U.0OD 


n 17c; 
u. 1 ( 


U.o4 It u.uo 


1 79 _|_ n 1 9 


2.40 ± 


0.24 


l.yo It U.4D 


0.670 


0.196 


0.49 ± 0.06 


1.68 ± 0.12 


1.54 ± 


0.24 


1.01 ± 0.42 


0.789 


0.218 


0.48 ± 0.12 


1.24 ± 0.20 


1.64 ± 


0.40 


0.92 ± 0.56 



4 Deconvolution Weighting 



Because of q 2 smearing, the sum of the P + weights in a given reconstructed 
q 2 bin will depend on both the underlying H+^q 2 ) for that bin as well as the 
underlying H+^^q 2 ) form factor for all other bins. 4 Hence the old Eq. 10 
correction becomes the non-local version given by Eq. 11 where for notational 
simplicity we just consider two q 2 bins. 

/ + (Y) = ip + . 13 = sSMl + i^SMl ? (11) 



where = ■ 1 M 1 is the sum of the Monte Carlo weights that reconstruct 
in q 2 bin 1 when generated in q 2 bin 1 and 1 w 2 = 1 P+ ■ 1 M 2 is the sum of the 
Monte Carlo weights that reconstruct in q 2 bin 1 when generated in q 2 bin 2. 
Eq. 11 can be generalized to Eq. 12. 



/+(V) 



V/+(V) 




2 



U>2 



\ 



( u2(Jl\\ 



HUqf 
\Hl{qD) 



(12) 



The solution to Eq. 12 is given by Eq. 13. 




v+^2 



W2 
+ Y2 



(H+)2 (ff+) 
2 wi 2 w 2 



-1 



l P + ■ X D 
2 P + ■ 2 D 



(13) 



4 Throughout this discussion we will use superscripts for reconstructed q 2 bin num- 
bers and subscripts for true q 2 bin numbers. 
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This solution is equivalent to weighting the data using the weights given by 
Eq. 14. 



-1 



TR2 



Mil W2 



1 1 
V 



2 P, 



(14) 



In order to obtain Hl(qf), for example, one weights each event by pi . By 
Eq. 14, pi is constructed from a sum of the vector 1 P + that contains the 
angular weights for events that reconstruct in q 2 bin 1 and the vector 2 P + 
that contains the angular weights for events that reconstruct in q 2 bin 2. 
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(b)i 
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i . . . . i . . . . i . . . 




q 2 GeV/c £ 



Fig. 5. Non-parametric form factors obtained for deconvolved data. The plots are: 
(a) H 2 + (q 2 ), (b) H 2 _(q 2 ), (c) H 2 (q 2 ), and (d) h (q 2 ) x H (q 2 ). 



Table 2 

Summary of deconvolved results. This is a tabular summary of the data of Fig. 5 
multiplied by 1000. 



q 2 


Hi 


H 2 


Hi 


h Q x H 


0.125 


0.508 ± 0.15 


1.24 ± 0.97 


17.00 ± 4.05 


4.72 ± 2.33 


0.375 


0.304 ± 0.14 


1.71 ± 0.35 


0.33 ± 1.17 


0.41 ± 0.95 


0.625 


0.270 ± 0.10 


1.73 ± 0.18 


2.82 ± 0.82 


2.37 ± 0.71 


0.875 


0.698 ± 0.12 


1.16 ± 0.16 


0.54 ± 0.49 


-0.14 ± 0.65 



Figure 5 and Table 2 show the result of a four q 2 bin deconvolution of the data. 
Only four bins are used since as the q 2 smearing exceeds the bin separation 
the "resolution" matrix that is inverted in Eq. 14 becomes increasingly more 
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singular resulting in greatly inflated error bars as well as strong negative corre- 
lations appearing between adjacent bins [5]. Essentially the same features that 
appear in the un-convolved representation of the data seen in Figure 4 appear 
in Figure 5 but on a coarser scale. There is a small overshoot in the first bin 
of both and Hq, but generally both the shape and relative normalization 
of the form factors are a reasonable match to the model of Reference [2] . 

We have estimated systematic errors by studying the stability of the results 
to changes in the number of angular bins, changes in the analysis cuts, and 
changes in our assumed background level. For the three form factor products 
describing the D + — > K component - H\(q 2 ), H 2 (q 2 ), and Hl(q 2 ) - the 

systematic errors are estimated to be less than 20% of the statistical error apart 
from the first q 2 bin for H^_(q 2 ), where we assess a systematic error equal to 
the statistical error. This first q 2 bin also requires a large 30%) background 
subtraction. For the h (q 2 ) x H (q 2 ) form factor product describing the s-wave 
interference we assess a systematic error equal to 30% of the statistic error. 
The background subtraction for the ho(q 2 ) x H (q 2 ) form factor product is 
also large (~ 50%) for the first three q 2 bins. The errors quoted in Tables 1 
and 2 do not include these systematic errors. 



5 Summary 

We presented the first non-parametric analysis of the helicity basis form fac- 
tors that control the decay D + — > K~ir + fi + is. We used a projective weighting 
technique that allows one to determine the helicity form factor products by 
weighted histograms rather than likelihood based fits. This method should 
prove to be a valuable technique for the e + e~ charm factories that have much 
better q 2 resolution. The non-parametric technique can also be used for other 
four body semileptonic decays like D + — > p and Df — > <fi £ + ue- We pre- 
sented both an un-convolved and deconvolved representation of the data. Both 
representations were a reasonable match to the spectroscopic pole dominance 
assumption with form factor ratios given in [2]. The q 2 dependence of the form 
factor governing the s-wave contribution to D + — > K^-K + ji + u is also studied 
for the first time. We find that the shape of the s-wave form factor (h (q 2 )) is 
reasonably consistent with the H (q 2 ) form factor as assumed in Reference [2]. 
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